X (3872) as a 1 D 2 charmonium state 
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The 1 I?2 charmonium assignment for the A" (3872) meson is considered, as prompted by a recent 
result from the BABAR Collaboration, favouring 2 quantum numbers for the X. It is shown that 
established properties of the A(3872) are in a drastic conflict with the 1 D2 cc assignment. 
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Seven years after its discovery, the A" (3872) meson has 
confirmed once again its reputation of enfante terrible of 
meson spectroscopy. The state was first seen by Belle [l[ 
in the tt + it~ J/ip mode and then confirmed, in the same 
discovery mode, by the CDF 0, D0 @], and BABAR Q 
Collaborations. According to the CDF analysis of the di- 
pion mass spectrum and the angular distribution in the 
TT + Tr~J/ip mode 5], only 1 ++ and 2 h assignments are 
able to describe the data. Then, while the nature of the 
state remains controversial, there are good phenomeno- 
logical reasons to assign it 1 ++ quantum numbers. 

First, the X resides at the D°D*° threshold, which 
prompts a considerable admixture of a molecule in its 
wave function. Furthermore, CDF concludes that the 
7r + 7r~ come from the p @ which, together with the Belle 
observation of the ojJ/ijj mode [7|, points to a consider- 
able isospin violation. The latter can be explained natu- 
rally in the molecular model of the A, which implies 1 ++ 
quantum numbers. In addition, the X was also observed 
in the D°D*° mode with a significant rate @-[loj]- Both 
pJ/ip and D°D* Q modes were analysed simultaneously 
in Refs. [llT - [l3j ]. and it was shown that indeed the data 
were compatible with a large admixture of the D°D*° 
molecular component in the wave function of the X. 

However, a recent analysis of the decay B — > KuiJ/ip 
data performed by the BABAR Collaboration [l4| indi- 
cates that inclusion of an extra unit of the orbital angular 
momentum in the ujj/ip system improves significantly 
the overall description of the observed 7r + 7r~7r° mass 
distribution, which implies a negative P-parity of the 
A (3872) state. Although this new BABAR result is fully 
compatible with the 2 h assignment allowed by CDF, if 
confirmed, it clearly challenges our understanding of the 
charmonium spectroscopy above the open-charm thresh- 
old. Here we investigate the most conventional explana- 
tion for the 2 h A (3872) as the l 1 ^ charmonium state. 

In case of charmonium D-levels we have an experimen- 
tal anchor at our disposal — the ■0(3770) vector state 
which is dominantly a cc state, with the angular momen- 
tum of the quark-antiquark pair L = 2 and the total 
quark spin S = 1 i}D 2 state has L = 2 and S = 0). 
As c-quark is heavy, the spin-orbit force, which splits 
spin-triplet and spin-singlet levels, is not large and all 
D-levels are degenerate in the leading-order approxima- 
tion. Hence, one may use the data on the 3 Di level to 
estimate the mass and matrix elements of the 1 D 2 level. 



It became clear long ago that the 1 D 2 assignment for 
the X disagreed with quark model mass estimates (see, 
for example, Refs. [R5l - [l7l ]). Indeed, quark models usually 
predict the 1 1 D 2 mass in the range 3770 ~ 3830 MeV, 
while the mass difference between the 1 1 D 2 and l 3 Di 
levels is predicted to be, averagely, 20 -7- 30 MeV. Thus 
quark models cannot accommodate the A (3872) as a 1 D 2 
state. The same conclusion was drawn in a recent paper 
[l8T ] . One might think that inclusion of various D-meson 
loops changes this statement. It is not the case, however. 
Loop calculations in the Cornell decay model [19j and in 
the 3 Pq decay model [2(| give for the mass of the 1 1 D 2 
level 3838 MeV and 3800 MeV, respectively. 

The arguments based solely on the mass calculations 
are, of course, not enough to rule out the charmonium as- 
signment for the X. However further reasons for the D 2 
interpretation to be problematic have started to show up. 
First, radiative decay transitions 1 D 2 — > 'yJIMib') rates 
are shown to be incompatible with the data [21j. Second, 
the production cross section of the 1 D 2 level at CDF is 
predicted to be much smaller than the one actually ob- 
served for the X [Ta j. In this paper we identify a couple 
of new problems. Namely, we expand on the issue of 
radiative decays and discuss the D°D*° mode of the X. 

The BABAR Collaboration has reported the following 
rates for the decays A(3872) -» 7 J/V>(V>'(3686)) 

Bi = Br(B ± ->■ K ± X)Br{X ->■ -yj/tp) 
= (2.8 ±0.8 ±0.2) x 10~ 6 , 

(1) 

B 2 = Br(B ± -> K ± X)Br{X -> 7?//) 
= (9.5 ±2.9 ±0.6) x 10~ 6 . 

In the meantime, the upper limit on the total branching 
fraction B ->• ATA (3872) imposed by BABAR [H is 

Btot = Br(B -> KX) < 3.2 x 10~ 4 . (2) 

Below we demonstrate that measurements (p} and ([2j) 
cannot be reconciled with each other under the assump- 
tion of the X being a 1 1 I?2 charmonium. To this end we 
notice that the leading multipole for the 1 D 2 — > 7V (V is 
a vector charmonium) transition is Ml, with the width 
given by a standard formula (see, for example, Ref. [17(): 

r( 2S+1 L. ; -> 2S ' +1 L'j^ (3) 
4 2J' + 1 ae 2 c . 2 3 

= 3 ar+r^'^^i^i^W 
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where m c is the charmed quark mass, e c = 2/3, ipi(ipf) 
is the initial(final)-state radial wave function, and E 1 
is the photon energy. In this formalism, the transition 
1 D 2 7^/0("0') i s a so-called hindered transition, so 
that (ipf\ipi) — sin 9, where 8 is the 3 Si — 3 Pi mixing 
angle. Thus the amplitude simply vanishes if J/ijj(ip') is 
assumed to be a pure 3 Si state. The standard value for 
the ip' is 9 ps 12°, which gives (for m c = 1.5 GeV): 

P ( 1 P 2 (3872) -> 7?//) « 6.6[keV] sin 2 8 ps 0.29 keV. (4) 

Notice that, being almost a pure 3 Si state, J/tp pos- 
sesses a tiny mixing angle 6, so that even a much larger 
photon energy (P 7 = 698 MeV for the jJ/tp final state 
versus E 1 — 186 MeV for 'yip') cannot provide a siz- 
able contribution of this, formally leading, Ml transi- 
tion. Therefore, contributions of higher multipoles have 
to be considered. In Ref. [2l[ both widths were calculated 
in a quite elaborated (though rather mo del- dependent) 
NRQCD approach, with the result: 

P ( 1 L> 2 (3872) -> 7 V/) ss 0.45 -j- 0.5 keV, (5) 
P ( 1 L> 2 (3872) -» 7 J/*0) w 6.8 — 9.5 keV, (6) 

claimed in Ref. HH to contradict the BABAR data {TJ. 

For the case of the -0', the Ml contribution from the 
3 Si — 3 Di mixing to the result ([5]) is indeed dominant 
and it is in a good agreement with the simple estimate 
(|4|. Notice that formula (j3j does not take into account 
recoil corrections, while the formalism of Ref. [2l[ ac- 
counts for the recoil only via the multipole expansion. 
Because of a small photon energy this seems reasonable 
for the ip' (3686) final state while, for the J/ ifi final state, 
the photon energy is much larger, so that the value (|6]) 
is probably an overestimation. 

In order to estimate the total width of the 1 D 2 charmo- 
nium we notice that, as is well-known (see, for example, 
Refs. [HI, Qj], [HI), the main radiative transition of the 
1 D 2 state is 1 D 2 -> 7^(3525) with the width: 

P ( 1 L> 2 (3872) -» 7 1 P 1 (3525)) ps 460(345) keV [H([2l|). 

(7) 

Alternatively, this width can be estimated from the 
measured branching fractions for the transitions [25| : 

Br (0(3770) -= 7 3 P (3415)) = (7.3 ± 0.9) x 10" 3 , 
Br (^(3770) -> 7 3 Pi(3510)) = (2.9 ± 0.6) x 10~ 3 . 

Indeed, the leading multipole is El, with the width: 
r ( 2S+1 Lj 2S ' +1 L' r ) = \c }i 8 SS >ela\ <^|r|^> \ 2 E 3 , 

where C fi = max(L, L')(2J' + { j' £ f } • In thc 

heavy-quark limit, the dipole matrix element (0>/l r l ri) is 
the same for all D — > P transitions. This gives 

r( 1 L> 2 (3872) 7 1 P 1 (3525)) ps 340 = 440 keV, (8) 



in a good agreement with Eq. 

Hadronic modes of the ^^P 2 charmonium were esti- 
mated in Ref. [H)]. These are the light hadrons ("<?<?") 
modes and the ?7 c 7r7r mode: 

P ( 1 £) 2 (3872) -> light hadrons) ps 190 keV, 
P ( 1 £) 2 (3872) -> ^tttt) w 210 ± 110 keV. 

Together with the radiative decay modes these give for 
the total width the value P tot ps 800 keV. In principle, 
r tot should also include a contribution of the DD* modes. 
However, as will be shown below, with the suppression 
of the D°D*° mode, this contribution is negligible. 

Therefore, if we use, in accordance with the Ref. [2l| 
calculation, the value of about 8 keV for the jJ/tp width, 
we get: 

Br{B -> KX) = (800/8) X B x = (2 3) X 10~ 4 , (9) 

which is compatible with Eq. ([2]l, In the meantime, a 
similar estimate for the decay X — > jip' gives: 

Br(B -> KX) = (800/0.5) x B 2 = 16 x 10" 3 , (10) 

where r(X — > jip') = 0.5 keV as per ([S]) was used. This 
value is awfully larger than the upper limit ([2]). To rec- 
oncile B 2 with Btot one needs to decrease r to t fifty times. 

We conclude in such a way that the data on the ra- 
diative decays of the X(3872) do not allow for its 1 D 2 
charmonium interpretation, if the BABAR result on the 
70' mode holds true. Notice, however, that the most 
recent Belle results [26| read: 

Br(B ± -> K ± X)Br{X -> jJ/ip) 

= (1.781^4 ±0.12) x 10" 6 , 
Br(B ± -> K ± X)Br(X -+ 7 0') < 3.4 x 10" 6 , 

which suggests that the BABAR and Belle measurements 
for the 70' mode contradict each other. 

Let us now consider the D°D°ir° mode. In 2006 
the Belle Collaboration reported an enhancement of 
the D D tt° signal observed in the reaction B + — > 
K + D°D°Tr just above the D°D*° threshold 0, at 
M x = 3875.2 ± 0.7±5ll ± 0.8 MeV, with the branching 

Br(B + -> K + D a D a ir a ) = (1.02 ± 0.3li°; 2 ^) x 10~ 4 . 

The peak was confirmed by the BABAR Collaboration 
as well 0. However, recently the Belle Collaboration 
announced a new analysis for the D*°D° case and a 
lower peak position was obtained than reported before, 
namely, M x = 3872.9±g;2±g;| MeV, with the branching 

Br(B + -> K + D°D*°) = (0.8 ± 0.2 ± 0.1) x 10" 4 . (11) 

This enhancement was associated with the V(3872) state 
seen in the D°D*° mode (here and in what follows an 
obvious shorthand notation D°D*° = D°D*° + D°D*° is 
used) , and the mass shift was attributed to the proximity 
to the D°D*° threshold. 
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To proceed we find the ratio of branching fractions: 

R = Br{B -> KD°D*°)/Br(B -> KX) > 0.25, (12) 

where the lower limit for R was deduced from the data 
quoted in Eqs. ([5} and (TTTI) . In what follows we argue 
that it is not possible to reproduce such a large value of 
the ratio R under the assumption of the X being the l 1 D 2 
charmonium. Indeed, it is claimed in Ref. [27J that the 
peak position in the DD* invariant mass depends on the 
orbital momentum / of the DD* pair. In particular, it is 
shown that with / = 1 it is quite easy to produce a peak at 
about 3 MeV above the DD* threshold, accommodating 



in such a way both BABAR [9] and old Belle [8J mea- 
surements. Depending on the model parameters, a peak 
much closer to the threshold can also be reproduced with 
I = 1, so that there is no contradiction with the new Belle 
data [l(| either. The value / = 1 corresponds to the 2 h 
quantum numbers of the X and therefore suggests the 
1 D 2 assignment for the latter. However, then the D°D*° 
rate behaves as k 3 (k being the relative momentum in the 
D°D*° system), so the proximity to the D°D*° thresh- 
old implies a considerable suppression of the production 
rate. Below we make this argument quantitative. 

The ratio R can be calculated as 



dM (dBr(B -> KD°D*°)/dM) 



/ dM (dBr{B -> Knon{D°D*°)/dM) + dM (dBr(B KD°D*°)/dM) 
JM- Jm- 



where the integration takes place over the mass region where the A(3872) resides, conveniently defined as M± — 
M ± 10 MeV, with M being the A(3872) mass. The D°D*° and non(D°D*°) rates entering expression JTSJ are: 

dBr(B -> KD°D*°) _ B g ( 1 £ , 2 -> D°D*°) k 3 dBr{B -> Knon(D°D* )) _ B r (non(L>° D* )) 

dM ~~ 27r (M — M ) 2 + TJji ' dM ~~ 2n (M — M ) 2 + lf o4 /4 ' ^ ' 



where Mth = m(D° D*°) , and the constant B absorbs the 
details of the short-ranged dynamics of the 6-quark decay. 
Due to the factor fc 3 the expression for the D°D*° does 
not take a Breit-Wigner form. To account for the finite 
width of the D*° we assume for k(M) a simple ansatz 

H3 k(M) = ^yW~M th ) 2 + r„ 2 /4 + (M - M th ), 

where \x is the D°D*° reduced mass, _T» ~ 65 keV is 
the width of the D*° meson estimated from the data 
[25| on the D*^ meson. The standard expression for the 
two-body relative momentum is readily reproduced as 
-T* — > 0. Finally, anticipating a strong suppression of the 
D°D*° mode, we substitute r(non(D°l)* )) w r tot . 

The coupling g ( 1 I?2 — > D°D*°^j can be estimated in 
the 1 Z?2 model for the X using the 3 Di state ^(3770) as 
a benchmark (pdd is the relative DD momentum and 
the charged-neutral meson mass difference is neglected) : 

r( 3 D 1 ^DD)=g( 3 D 1 -+DD)p i DD . (15) 

We now invoke the "loop theorems" proven in Ref. [IK . 
In particular, it is shown in this paper that, in the 
heavy-quark limit, strong open- flavour total widths for 
the states in a given {NL} multiplet (N is the radial 
quantum number while L is the quark-antiquark orbital 
angular momentum) are equal. The heavy-quark limit 
implies that (i) the initial states are degenerate in mass 
and have the same wave functions within a given multi- 
plet and (ii) the final two-meson states exhibit the same 



degeneracy. The decay model should satisfy some general 
conditions listed in Ref. [29| (for example, the popular 
3 Po pair creation model satisfies these conditions, and so 
does the Cornell decay model). 

Specifically, in the ideal heavy-quark world, the masses 
of all ID states are identical, and the masses of the final- 
state D and D* mesons are identical too. The partial 
widths into certain D^D^ channels depend on quan- 
tum numbers of a given initial state, while the sum of 
partial widths over all possible D^*'D^*' final states is 
the same within a given ID multiplet. In the real world, 
if the quark-antiquark pair in the initial meson is heavy, 
the theorem is violated mainly by spin-dependent inter- 
actions, which remove the mass degeneracy both in the 
initial and final states. One may write therefore: 

g( 1 D 2 ^D°D*°)=g \C( 1 D 2 )\\ 

g( 3 D 1 ^DD)=g \C( 3 D 1 )\ 2 , 

where go is the coupling constant common for all mem- 
bers of the ID multiplet, while C ( 1 D 2 ) and C ( 3 L>i) 
are the spin-orbit recoupling coefficients for the 1 D 2 — > 
D°D*° and 3 D\ -» DD decays, respectively (notice that 
both charged and neutral DD channels contribute to the 
coefficient C ( 3 -Di), while only the D°D*° channel con- 
tributes to the coefficient C{^D 2 ^). These spin-orbit 
recoupling coefficients were calculated in the Cornell de- 
cay model (see Table II of Ref. [19() and in the 3 Pq decay 
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TABLE I: The ratio R (see Eq. HU)) for various values of the 
X(3872) mass Mo and the total width r to t- 



M , MeV 


r tot = 200 keV 


r tot = 800 keV 


r tot = 3200 keV 


3870.8 


0.023 


0.023 


0.022 


3871.4 


0.033 


0.032 


0.029 


3872.0 


0.073 


0.052 


0.038 



model (see Table IV of Ref. Both models yield 

\C( 1 D 2 )\ 2 = | |C ( 3 i^i) | 2 , so that, since the DD mode 
is dominant for the "0(3770), we estimate the coupling 
g ( x Da -> D°D*°) as: 

g^^jfiD^^W^m. ( 16) 

In Table |H we list the results for the ratio R for 
several values of r tot . The masses used are m(D°) = 
1864.84 MeV and m(D*°) = 2006.96 MeV. The V(3770) 
width is 23 MeV. Finally, for the mass Mg we take the 
same values as used in Ref. [27j. As described above, 
r tot = 800 keV is our preferred value. We have also cal- 
culated the ratio R for r tot four times smaller as well 
as four times larger than 800 keV, the latter value be- 
ing a bit larger than 2.3 MeV quoted in PDG [25[ as the 
upper limit for the width of the X. Clearly all values 



of R listed in the Table Q] are far too low in comparison 
with the value (|12|) deduced from the data. In addition 
we confirm the D°D*° lineshapes obtained in Ref. [27| . 
however, the rate appears to be quite small. Thus we 
conclude that the data on the D°D°ir° mode contradict 
the D2 charmonium interpretation of the X. 

To summarize, we have shown that the 1 1 £*2 char- 
monium assignment for the A(3872) meson contradicts 
the existing data on its radiative decays and its D°D n° 
mode. Our study does not challenge the 2 ^ quantum 
numbers. We rather claim that, if the aforementioned 
experimental data are taken as a true guide, the conven- 
tional charmonium model is not able to accommodate 
for the 2"+ AT(3872). If the BABAR result on the quan- 
tum numbers of the X(3872) persists, it would mean that 
some kind of a new interloper enters the game. 
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